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ABSTRACT 

In quantum theory of gravity, we expect the Lorentz Invariance Violation (LIV) and the modification 
of the dispersion relation between energy and momentum for photons. The effect of the energy- 
dependent velocity due to the modified dispersion relation for photons was studied in the standard 
cosmological context by using a sample of Gamma Ray Bursts (GRBs). In this paper we mainly 
discuss the possible LIV effect by using different cosmological models for the accelerating universe. 

Due to the degeneracies among model parameters, the GRBs’ time delay data are combined with the 
cosmic microwave background data from the Planck first year release, the baryon acoustic oscillation 
data at six different redshifts, as well as Union2 type la supernovae data, to constrain both the model 
parameters and the LIV effect. We find no evidence of LIV. 

Keywords: gravitation— (stars:) gamma-ray burst: general— (cosmology:) dark energy 


1. INTRODUCTION 

Above the Planck energy scale Epi, we expect a quan¬ 
tum theory of gravity in place of Einstein’s general rel¬ 
ativity. The quantization of space-time will lead to the 
modification of the dispersion relation between energy 
and momentum of a particle with mass m and the break¬ 
down of Lorentz invariance. In some quantum gravity 
models, the consequence of the Lorentz Invariance Vi¬ 
olation (LIV) is t he energy dependence of the speed o f 
light in vacuum (|Amelino-Gamelia et al. I B9981 I200II) . 
In these scenarios, the energy dependent velocity of light 
is u = c(l — ^E/Eqg) with effective quantum gravity 
energy scale Eqg, the speed of high energy photons is 
slower, and low energy photons will reach us earlier. 
Therefore, the measurement of the light speed in vac¬ 
uum can be used to study the LIV effect. 

In the past years, both astrophysics and particle 
physics’ experiments have been proposed to test the 
LIV and quaiitum gravity effects dM attingl v et al. II2005I : 
iSarkar et al. 1 1200^ lAmelino-Gamelia. II20R{[) . ^cause 
photons with different energy will reach us at dif¬ 
ferent time, and the energy of Gamma-Ray Bursts 
(GRBs) which are the most luminous explosions in 
the universe with short duration ranges from KeV 
to GeV, GRBs with intense bursts of 7 -ray photons 
from cosmological distances are proposed to measure 
the difference of the arrival time of photon s with 
different energy (lAmelino-Cam elia et al.~1 B998L 120021 : 
iKowalski-Glikman &: Nowak 1 1200211. GRB s were also 
used to constrain LIV theories (lEllis et al. 112001^120061 : 

Rodriguez Martinez fc Piran 1120001 . 1.Iacob et al. 1120071 1. 

Biller et al. I ( 1999f l used the data for a TeV 7 -ray flare 
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coming from the active galaxy Markarian 421 to set a 
bound on the energy scale of quantum gravity as Eqq > 

6 X IQi® GeV. lAbdo et al. I (I2nn9ell used GRB 090510 
to study the LIV effect and get Eqg > 1.2i5pi, the en¬ 
ergy scale of quantum gravity was found to be Eqg > 

1.3 X IQi® GeV by using GRB 080916C. The LIV effects 
due to extra dimensi ons can also be constrained by GRBs 
l|Baukh et al. |[200^ . Experiments with cold atoms were 
also proposed to constrain the modified dispersion re¬ 
lation due to quantum gravity (jAmelino-Gamelia et al.l 

Hool). 

In order to measure the LIV, statistical and possible 
systemat ic uncertainties must b e minimized. For this 
purpose. lEllis et aT~l (|2003l . l200^ developed a method to 
analyze samples of GRBs with different redshifts and en¬ 
ergy bands. The technique has the advantage of extract¬ 
ing time-dependent fe atures from the signals of many 
GRBs. In the analysis (| Ellis et al. Il2003fl . they used both 
the BATSE and OSSE data from the Gompton Gamma 
Ray Observatory 0. By adding larger samples of GRBs 
with known redshifts from HETE 0 and SWIFT 0 , the 
observed time delay was formulated in terms of linear 
regression where the intercept denotes intrin sic time de¬ 
lay and the linear term d enotes LIV effect (|Ellis et al. I 
120081 : I Jacob et al. 1 1200811 . By using the concordance 
AGDM model with Ba = 0.7, they found no strong evi¬ 
dence of LI V, and they obtained Eqg > 1.4 x 10*^® GeV 
(|E11is et a,l. 112001 . 

Since the distances from the sources to the observer 
depend on cosmological models, the conclusions on the 
LIV and quantum gravity effect may be affected by the 
cosmological model used in the analysis. However, only 
the concorda nce AGDM scenar io {0\ = 0.7, Dm = 0.3) 
was studied (lEllis et al. Il2003l l. Therefore, it is neces- 
sary to conside r other cosmological models. Recently 
iBiesiada et al. I (12009( 1 studied the LIV by using the 
quintessence and Chaplygin Gas model and they found a 

^ ftp://legacy.gsfc.nasa.gov/compton/data/batse/ascii/data/64ms/ 

® http://space.mit.edu/HETE/Bursts/ 
http://swift.gsfc.nasa.gov/docs/swift/archive/ 
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weak evidence of LIV. In their analysis, the degeneracies 
among model parameters were neglected, as they were 
taken as fixed values. In this paper we mainly discuss 
the possible LIV by using different cosmological mod¬ 
els which explain the present cosmic acceleration. Com¬ 
pared with the previous work, in addition to the rele¬ 
vant parameters quantifying LIV, the cosmological pa¬ 
rameters are also treated as free parameters. We adopt 
the Markov Chain Monte Carlo (MCMC) technique to 
constrain LIV and model parameters with the observa¬ 
tional data. In order to derive tighter constraints on the 
model parameters, we combine the time delay data from 
GRBs with the cosmic microwave backgro und (CMB) 
data from the Planck first year release (|Ade et al. I 
l20l4rWang et al. 1120131). the bar y on acoustic osci l lation 
(BAO) data llBlake et al. I l201ll: iBeutler et al. I 120111 : 
iCong et al. I I20l4 Percival et al. 1 1201011. and the 557 
Union2 SNela data ( Amamillah et al. II 201 I . This paper 
is organized as follows: we introduce the LIV in section 
11. The GRBs’ time delay and other observational data 
are discussed in section III. The constraint results on LIV 
and cosmological parameters with different cosmological 
models are presented in section IV, and finally, the con¬ 
clusions are drawn in section V. 

2. THE LORENTZ INVARIANCE VIOLATION 


The modified dispersion relation due to qu antum grav¬ 
ity models can be paramet erized as follows (|Ellis et al. I 
I 2 OO .4 iBiesiada et al. |[2009[ 1 


E'^ -p^c^ = eE 


E y 

inEga J 


( 1 ) 


where is a dimen sionless parameter, = 1, ^2 = 10“^ 
(| Jacob et al. 11200^ . e = -|-1, and the effective energy 
scale Eqc of quantum gravity is expected to be near 
the Planck scale. Because the LIV e ffects are small, 
(|Ellis et al. 1120081 IBiesiada et al. lU009ll . we consider the 
n = I term only in this study. 

Since t he arrival time for photons with energy E is 
equal to (|Jacob et al. 11200^ 120081 : IBiesiada et al. 1120071 

[2?10l 


thiv 




E 

Eqg 


(1 + ^') 


dz' 


( 2 ) 


so the the arrival times of two photons with different 
energy will be different, and the difference in time for 
the energy difference AA is 




AE r {l + z')dz' 
HqEqg Jo Hz') 


(3) 


where the dimensionless Hubble parameter Hz) = 
H{z)/Ho, Hq is the Hubble constant and H{z) is the 
Hubble parameter at redshift z. In order to account 
for the unknown intrinsic time lags, we fit the measured 
time-lags by including a parameter b specified in the rest 
frame of the source. The arrival time delays are fitted b y 
the formula Atobs = AIliv + b{l + z) (|Elhs et al. II2006I) . 
Therefore, the simple linear fitting function is 


—2^ = a^ivK + b, (4) 


where the intercept represents the intrinsic time lags, 
the LIV effects are encoded in the slope aLiv = 


AE/(Hq Eog) which is related to the scale of Lo rentz 
violation (I Jacob et al. Il200§: IBiesiada et al. 1 1200711 and 
K = {I + z)~^ fj dz'(l + z')lk(z') is related to the mea¬ 
surements of cosmic distances. The model dependence is 
through the function K(z) which will be calculated for 
three popular cosmological models. 


3. THE OBSERVATIONAL DATA 

In this paper we use the BATSE data with 9 light 
curves whose time resolution is 64 ms and redshifts span 
from from z = 0.835 to z = 3.9, the HETE data with 15 
light curves whose time resolution is 164 ms and redshifts 
span from z = 0.168 to z = 3.372, and the SWIET data 
with 11 light curves whose time resolution is 64 ms and 
redshifts span from z = 0.258 to z = 6.2 9, the data 
was shown in Table 1 of lEllis et al. I (j2006[ l. So we use 
the data of time lags between different energy channels 
measured from the light curves of 35 GRBs with re dshifts 
from z = 0.168 to z = 6.29 (|Elhs et al. 1 120061) . The 
spectral time lags are obtained from the light curves in 
the 115-320 keV energy band with respect to those in the 
lowest 25-55 keV energy band. To test the LIV effect, we 
apply the statistics, here is 


Xgrbs 


NgRBs 

E 


2=1 


Ati Atobs 
o-i 


(5) 


where Ati and Atobs respectively denote the theoretical 
and observational values of time delays of GRBs, and <7^ 
is the observational uncertainty. 

In fitting the data, in addition to the slope a^jy and 
the intercept 5, we also need to fit the model parameters. 
In previous studies, the model parameters are fixed and 
the degeneracies among parameters are neglected. To 
account for the degeneracies among these parameters, it 
would be better to treat the model parameters as free 
parameters and estimate their nominal values from the 
observational data. With this aim, we combine the CMB 
data from the Planck first year release, the 557 Union2 
SNela data, and the BAO data with the data of spectral 
time lags from GRBs to constrain the model parameters. 
Since we have at least four p arameters to fit, we take the 
MCMC (iLewis et al. 1121)?)^ technique to constrain the 
model parameters. 

For the 557 Union2 SNela data, the distance modulus 
Hz) is measured at different redshifts, and the theoreti¬ 
cal value of the distance modulus is 

/r = 51og-^-k 25 = 51ogioi7odL - Mo, (6) 

Mpc 

where Mo = 5 log^Q [Ho/ (100 km/s/Mpc)] — 42.38, and the 
luminosity distance is c?l = (1 -k z)Hy^ fj dz'/E{z'). To 
fit the 557 Union2 SNela data, we calculate 


2 _ ^ [Hz^) - ^J'obsiz^)y 

XSNe - 2^ -2 

i=l 


(7) 


where Hzi) and Pobsizi) are the theoretical and observed 
distance modulus for the SNela at redshift z^. The cr^i is 
observational error. For the nuisance parameter Hq, we 
marginalize it with a flat prior. 

The BAO A data consist of the measurements of A at 
three redshifts z = 0.44, z = 0.6 and z = 0.73 from Wig- 
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Table 1 

The corresponding inverse covariance matrix of A. 


Redshift slice 

0.2 < 2 < 0.6 

0.4 < 2 < 0.8 

0.6 < 2 < 1.0 

0.2 < 2 < 0.6 

1040.3 

-807.5 

336.8 

0.4 < 2 < 0.8 


3720.3 

-1551.9 

0.6 < 2 < 1.0 



2914.9 


gleZ dark energy survey (|Blake et al. 1120111 : iGong et al. I 
I2014f) . and the parameter A is defined as 


^ (s) 

Z 

wher e the effective distance is given by (lEisenstein et al. I 


Dv{z) = 




-I 1/3 


(1 + 2)2 E{z) 


(9) 


The value for this data set is 

Xbaoi = ^AbaqC^bao ^AAbaO) (10) 


where Cbao”^ is the corresponding inverse covariance 
matrix (See Table [T|). 

In addition to the above A parameter, we also con¬ 
sider the BAO distance ratio dz at 2 = 0.2 and 2 = 
0.35 from SPSS dat a release 7 (DR7) galaxy sample 
(|Percival et al. 11201^ a nddn.infi = 0.336+0. 015 from the 
6 dFGS measurements (iBeutler et al. l[20Tl[l . The BAO 
distance ratio is 


, _ T,{za) 


( 11 ) 


where the comoving sound horizon is 


Ts 



Ca{z)dz 

E{z) 


( 12 ) 


the sound speed Cs{z) = l/-\/3[l + i?h/(l + 2 )], and 
Rb = 30;,/i2/(4 X 2.469 x 10~^). As usual, we fi t the 
drag redshift Zd as follows (|Eisenstein et al. N1998I1 


Zd = 



1291(O„/i2)0-25i 
1 +0.659(O^h2)0-828 

0.313(G„h2)-o-4i9[l - 

0.238(11™/i2)O-223^ 


[l + 6 i(fi!h/i 2 )'^], 

b0.607(O„/l2)0'674]^ 


(13) 


The value o f the BAO distance ratio is 

(|Percival et al. 11201011 

Xbao2 = ^Pbao^bao ^APbao, (14) 

where APbao = Pth — Pobs, Pobs is the observed dis¬ 
tance ratio, and Cbao is the covariance matrix for the 
distance ratio. To use the BAO data, we minimize 


2 _ 2 ,2 , ( 1 I 0.106 — 0.336)2 

Xbao — Xbaoi + Xbao2 + ' 


0.0152 


(15) 


For the CMB measurement, we use the derived pa¬ 
rameters including the acoustic scale Ig, the shift param - 
eter R, and VLbh? (|Ade et al. 112014 iWang et al. 11201^ 
from the Planck temperature measurements combined 


Table 2 

The Icr and 2cr constraints on a and h for different cosmological 
model. 


Cosmological model 

Regression coefficient a 
and intercept b 

ACDM 

_ Q „+ U . U ^ i /,+ U . 141 b 

a— 0718,-0.1415 

_0 qqqio+0.0154,+0.0308 

0 — U.UUUiC)o 0155.-0.0305 

u;CDM 

^ ~ 0702,-0.1392 

u r. nnni c+0.0153,+ 0.0303 

b = -0.00015 o o, 54+0 0304 

CPL 

^ _ n r»-| QO + U-U fii-ii + U . 14 

a — u.uioo_0.0711,-0.1401 

h n nnni o+0.0159,+0.0306 
b — U.UUUi«_o.0155.-0.0301 


with leasing, as well as WMAP polarization data at low 
multipoles with / < 23. The acoustic scale is 

, St dz/E{z)] 

= ^ -7771-’ 

where rs( 2 *) = E[o~^ Cs{z)/E{z)dz is the comoving 
sound horizon at the recombination. The shift parameter 
is 


R{z:Y = ^{^ sinn (f . (17) 

y/Wk\ \ Jo E{z)J 

The recombination redshift 2 * is fitted by (|Hu et al. I 
[19^ 

2 * = 1048[1 + 0.00124(nb/i2)-o™][l + 

0.0783(06h2)-o-238 
l + 39.5(ll6/l2)0.763 ’ 

0-560 

l + 21.1(llbh2)l-81- ^ ^ 

The x^ value for the CMB data is 

XCMB = ^PcMB^CMB ^APcMB, (19) 

where APcmb = Pth - Pobs, Pobs and Pth are the ob¬ 
served and the theoretical values of the derived parame¬ 
ters, respectively, and Cumb is the covarian ce matrix for 
the derived CMB data (IWang et al. 11201^ . 


4. CONSTRAINTS ON THE LIV PARAMETERS AND 
ANALYSIS 

To explain the cosmic acceleration, an exotic energy 
component calle d dark energy with negative pressure 


19991: lAstier et al. 1 120061: Hicken et al. 1 

2009; 

Amanullah et al. 1 20101: iSnergel et al. 20031 

2007; 

Komatsu et al. 120091 120111; Tegmark et al. 1 

2004; 

Eisenstein et al. 1 120051: Cao et al. 

[Ml [2m 

L2allb: 

Gao & Gong 112014 IGong et al. 112013 

, |2014[). The most 


simple candidate for dark energy is the vacuum energy 
known as the cosmological constant A. Because of 
the huge difference between predicted and measured 
values of the vacuum energy, many other dynamical 
dark energy model s have also been considered, includ¬ 
ing quintessence (iRatra et al. I 1988 : iGaldwell et al. I 
1199811. phantom (iCaldwell et al. 1 12002,1200311. k - essenc e 


( Armendariz-Picon et al. I 120(31 


Chiba et al. I 1200211 


as well as quint om models (|Feng et al. I l200d 120061 : 
I Guo et ain 1200511 . In this paper, we consider three 








































































































4 


Pan et al. 



Figure 1. The marginalized Icr and 2a contours and the distri¬ 
bution of the LIV parameters for the ACDM model. The red cross 
represents the best-fit point. 

different dark energy models: the ACDM model, the 
dark energy model with constant equation of state 
parame ter w, and the Chevall i er-Polarski-Linder ( CPL) 
model (jChevallier et al. II2001I : iLinder et al. 1120031) . 

For the ACDM model, the equation of state parameter 
w = p/p = —1 and the Friedmann equation is 


h^{z) — + Oa, 


( 20 ) 


where is the matter energy density, and the dark 
energy density is Da. Since Dm + Da = 1 for a flat 
ACDM model, so we have only one independent param¬ 
eter Dm in this model. By fitting the ACDM model 
to the above data, we get the la and 2tT constraints: 

_ pv pv-I «-|-0.0717,-|-0.1416 7 p\ p\p\p\-j qTO.0154,-|-0.0308 

® ~ ‘-0.0718,-0.1415J " ~ 0155,-0.0305’ 

o n onT0.010,-|-0.022 i jj- r-l-0.9,-1-1.8 't'i,,.., ^ 

— ^*^^—0.011,-0.022’ o,nci siQ — 69.5_g g 7 ' -Lnc rc- 
suits show that Lorentz invariance is consistent with the 
data. We show the results on a and b in Table m The 
ID probability distribution of each parameter and the 
2D confidence contours for the parameters are shown in 
Fig.m 

For the flat icCDM model, the Friedmann equation is 


h^iz) — Dm(l + z)^ -l- Da:(1 + 


( 21 ) 


where the dark energy density Dx = 1 — Dm. Using the 
MCMC method, the marginalized Icr and 2a constraints 
on the model parameters are: a = —0.0168lg Q 7 Q 2 ’lQ i 3 g 2 , 

L _ ri nnni r+0.0153,-1-0.0303 _ pv qqo+0-011,+0.022 

o — —_Q Q 3 Q 4 , \Lm — - 0 . 021 ’ 

1 nc^+0-04,+0.079 i tt- r)+l.l,+2.1 ^ 

w — ~1.05_q Q 45 _Q Qgg, and Hq — 70.2_4 ^ _2 4 . There 
is no evidence of LIV in wCDM model. The results are 
shown in Table [2l and the marginalized plots are shown 
in Fig. [2j 

For the CPL model, the equation of state parameter 
is parameterized as w{z) = wq + Wiz/{1 + z) and the 
Friedmann equation is 

h^iz) = Dm(l+z)+Dx(l+z)3(^+“'“+“^) exp • 

m) 

In this model. Dm, wq, and Wi are the model pa¬ 
rameters. Fitting the CPL model to the combined 
data, we obtain the marginalized Itr and 2 ct constraints 
on the model parameters: a = —O.OISSIq q^^^’Iq+q^, 

L _ p\ p\p\p^i q+0.0159,+0.0306 q _ pA c^oo+0.012,+0.024 

U — '9.UUU10 _q Q 255 _o.0301’ 012,-0.022’ 


Figure 2. The marginalized la and 2a contours and the distri¬ 
bution of the LIV parameters for the uiCDM model. The red cross 
represents the best-fit point. 




Figure 3. The marginalized la and 2cr contours and the distri¬ 
bution of the LIV parameters for the CPL model. The red cross 
represents the best-fit point. 


-I pvr)+0.123,+0.246 +0.593, +1.053 i 

Wo — ~l-02_o 122,-0.237’ “ ~Q*^03 _q 595 _i 281’ 

Ho = 70.3l);o^I_+57®. We see no evidence of LIV in the 


CPL model. The constraints on a and b are shown in 
Table [21 the marginalized probability distributions and 
the la and 2cr contours of a and b are shown in Fig. [31 
In the previous analysis, the model parameters take 
fixed values and a weak indication of LIV was found 
(j Ellis et al. 1 120061: iBiesiada et al. 1 120091) . Here we con¬ 
sider the possible degeneracies among all the parame¬ 
ters and take all the parameters as free parameters, we 
find no evidence of LIV and all the results are consistent 
with each other. The conclusion is independent of the 
background model. Although the uncertainties on the 
parameters a and b become larger due to more fitting 
parameters, the difference on a for different cosmological 
model is still relatively small. To better understand the 
effect of the background model, we show the re-scaled 
spectral time lags Atobs/{^ + z) versus the K{z) function 
for the three different cosmological models in Figs. |4]|6l 
the data from GRBs is also shown in the figures. Because 
the background cosmological model changes the value of 
K{z), so the value of a also changes . Unlike the models 
considered in IBiesiada et al. 1 (|2009D . the wCDM model 
and CPL model considered here is close to the ACDM 
model as seen from the above best-htting values of w, wq 
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Figure 4. The re-scaled spectral time-lags Ato6s/(l + ^') versus 
K{z) for the ACDM model with Qrn = 0.29. The solid line repre¬ 
sents the best fit with the linear regression. 



Figure 5. The re-scaled spectral time-lags Aio 5 s/(l -h z) versus 
K{z) for the itjCDM model with flrn = 0.288 and w = —1.05. The 
solid line represents the best fit with the linear regression. 



Figure 6. The re-scaled spectral time-lags Atobs/{^ + z) versus 
K{z) for the CPL model with Qrn = 0.288, mo = —1.02 and mi = 
-0.203. 


and wi^ so our results are all consist with each other. 

The redshits of the GRBs span from z = 0.168 to 
z = 6.29, so we also consider the effect of the redshift 


Figure 7. The constraint on the slope a for the mCDM and CPL 
models normalized to ACDM model by using four sub-samples of 
GRBs. The GRBs are divided into 4 groups with upper boundaries 
at 2 = 1.0, z = 2.0, 2 = 3.0 and 2 = 6.3. The circles represent 
the CPL model and the triangles are for mCDM model. The red, 
black, blue and green are for the fits to data with redshifts 2 < 1, 
2 < 2, 2 < 3 and 2 > 3, respectively. 

distribution on the results. We divide the GRBs into 
4 groups with upper boundaries at z = 1.0, z = 2.0, 
z = 3.0 and z = 6.3. The first group has 11 GRBs with 
redshifts 0 < z < 1.0, the second group has 21 GRBs 
with redshifts 0 < z < 2.0, the third group has 27 GRBs 
with redshifts z < 3.0 and the fourth group contains all 
35 GRBs. We then fit the cosmological models to each 
group of GRBs, and the constraints on the slope a nor¬ 
malized to the ACDM model are shown in Fig. [T] The 
difference between different models are negligible in the 
first and the third redshift groups. For both wCDM and 
the CPL models, the effects of redshift distributions are 
similar although the deviation is a little larger for the 
CPL model. The biggest contribution comes from the 
GRBs with redshifts in the range 1 < z < 2 and z > 3. 
As we discussed above, the CPL model deviates more 
from ACDM than wCDM model does, so the value of 
the slope a for the ACDM model is more negative. The 
reason of the redshift-dependence is due to distributions 
of the GRBs. In the first and third groups, the average of 
the spectral time lags is close to zero. But we have more 
negative At data in the redshift intervals 1 < z < 2 and 
z > 3, therefore we see bigger deviations for the data up 
to redshift 2 and for the whole data. 


5. CONCLUSIONS 


We test the LIV with the 35 GRBS, the Union2 SNela, 
the CMB and BAG data for three different cosmolog¬ 
ical models. The slope a and the intercept b in the 
linear regression method as well as the model parame¬ 
ters are fitted to the combined data. For the ACDM 


model, the marginalized Icr and 2a results are: a = 

A A1'7~I”9.0717,-|-0.1416 ^ j j, a aaai 0154,-1-0.0308 

— 0.017_g _Q 2415, and b — — 0.00013_g Qigg _q 0305. 

For the wCDM model, the marginalized la and 2a 

results are: a = —O.OlGS^g'gygVgiQnn, and b = 

n nnni .+0.0153.+0.0303 ^ 1 

—0.00015_g gi 54 ’_o,g 3 g 4 . Lor the GPL model, we 


obtain the marginalized la and 2a results a = 

A A1 oo”I”6'9712,-1-0.14 ,, J l a aaai 0”I”6'9159,-1-0.0306 

— 0.0183_q 0711 _Q 1401, and b — O.OOO18_0 gigg _q gggi. 


These results are also summarized in Table [H Because 
the slope a is consistent with a = 0 for all the three mod- 
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els, there is no evidence of LIV. The results for all the 
three models are also consisten t. Our conclusi o n is i n 
conflict with previous result bv iBiesiada et al.~l (I2009D . 
In their studies, the model parameters are fixed at cer¬ 
tain values, so the degeneracies among parameters are 
neglected. The dynamical dark energy models with the 
fixed parameter values discussed in the previous work are 
very different from ACDM model, so they found different 
results for different cosmological models. In our analysis, 
we treat the model parameters as free parameters, and 
the model parameters are fitted with the latest SNela 
and BAO data. The best fitted wCDM and CPL models 
are close to ACDM model. Although the uncertainties 
on the model parameters becom e almost double com- 
pared with the previous analysis (|Biesiada et al. 11200911 
due to more fitting parameters, the differences between 
different model are relatively small and all the results 
are consistent. Because the GRBs’ data is not uniformly 
distributed, especially in the redshift intervals 1 < z < 2 
and z > 3, we see the small differences between differ¬ 
ent models come from the GRBs at the redshift intervals 
I < z < 2 and z > 3. In conclusion, our results show no 
evidence of LIV in the three cosmological models. 

The results are also in good agreement with the re¬ 
cent neutrino analysis, which discussed LIV using two 
cascade neutrino events wit h energies around I Pe V re¬ 
cently detected by IceCube iBorriello et al. I (I2013D . Yet 
it is worth noting that experimental probes of LIV are 
limited by the scarcity of GRB data. Other high-energy 
astrophysics experiments such as the pho ton time-delay 
meas urements from objects like Pulsars (lKaaret_et_jI^ 
1999 1 and Active Galactic Nuclei(AGN) (|Albert et al. 
200811 may provide complementary probe of LIV effect. 
Further studies are still needed to draw a more quanti¬ 
tative conclusion. 
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